* factor analysis replication

use MeasureComparisons.dta, clear
 
pwcorr measure5 polity2 aclpdemoc inverse_mean, sig o
polychoric measure5 polity2 aclpdemoc inverse_mean
matrix r = r(R)
global N = r(sum_w)
set more off
factormat r, n($N) ml protect(500)
